Slow dynamics of Zero Range Process in the Framework of Traps Model 
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The relaxation dynamics of zero range process (ZRP) has always been an interesting problem. In this study, 
we set up the relationship between ZRP and traps model, and investigate the slow dynamics of ZRP in the 
framework of traps model. Through statistical quantities such as the average rest time, the particle distribution, 
the two-time correlation function and the average escape time, we find that the particle interaction, especially 
the resulted condensation, can significantly influence the dynamics. In the stationary state, both the average 
rest time and the average escape time caused by the attraction among particles are obtained analytically. In the 
transient state, a hierarchical nature of the aging dynamics is revealed by both simulations and scaling analysis. 
Moreover, by comparing the particle diffusion in both the transient state and the stationary state, we find that 
the closer ZRP systems approach the stationary state, the more slowly particles diffuse. 
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The investigation of diffusion process plays a very impor- 
tant role in exploring the structures of systems and uncovering 
the physical mechanisms of dynamics in particle systems [ 1- 
4]. Recently, with the booming research results in the field 
of complex networks, particle diffusion in complex networks 
has once again become a hot problem 1I5I- I20I1 . For example, 
Noh et al. studied random walks in scale-free (SF) networks 
and found that the ratio between a coordination number and a 
characteristic relaxation time for each node essentially deter- 
mines the MFPT Hill. 

A common feature of above studies is that there is no in- 
teraction among particles at each node. Actually, particle in- 
teraction is unavoidable and can be found in many real net- 
works where it plays an important role in the dynamical pro- 
cesses in networks. One simple way to introduce the in- 
teractions among stochastic particles is so-called zero range 
process (ZRP), which has been recently adopted to investi- 
gate particle condensation in complex networks Il2ll - l25ll . It is 
shown that this condensation transition appears in a number of 
unexpected contexts such as wealth condensation in macroe- 
conomics l26ll . jamming in traffic ll27l - [30ll . coalescence in 
granular systems 1I3TI 13211 . and gelation in networks ll33l l34ll . 
To understand the influence of ZRP interaction on the par- 
ticle diffusion, we made an important first step to study an- 
other aspect of ZRP in scale-free networks, i.e., the diffusion 
features ll35ll . We found that the statistical quantities of dif- 
fusion can be significantly reduced by the condensation and 
can be figured out by the rest time of a particle staying at a 
node. In addition to these statistical quantities in the station- 
ary state, however, the relaxation dynamics of ZRP are fully 
overlooked. 

Recently, a traps model with interaction, in which parti- 
cles are attracted at each node (i.e. local minima) with the 
potential-energy landscape HU, was put forward to study the 
connection between the network of the potential-energy land- 
scape and the glassy dynamics 13711 . At low temperatures, the 
relation between the energy and the degree of a minimum can 



result in the slow dynamics of glassy systems. This provided 
a systematic integration of tools and concepts to investigate 
how network structures impact the particle diffusion. In this 
letter, we will try to set up the relationship between ZRP and 
traps model, and then investigate the slow dynamics in the 
ZRP interaction in the framework of traps model. 

Firstly, we briefly introduce the ZRP model. In complex 
networks, ZRP interaction means that a particle interacts only 
with other particles staying at the same node. Suppose N par- 
ticles are randomly put in a network with L nodes and each 
node i can be occupied by any integer number of particles 
n, = 0, 1, 2, • • • , N. Due to interaction, usually only part of 
the particles at a node can jump out and hop into its neighbor- 
ing nodes, which makes the particles be redistributed among 
all nodes. We divide the diffusion process of particles into 
two steps. In the first step, some particles at node i will jump 
out because of the interaction among the particles. Suppose 
a particle at node i jumps out at the rate p(rii) = n\ with 
8 € [0, 1] I2H - I23I1 . where the parameter S can be used to char- 
acterize the nature of interaction among particles at the node. 
In the second step, each jumping particle hops from the node 
i to one of its neighbors j randomly. 

In the stationary state, there is a condensation threshold 
5 C = 1/(7 — 1), and a finite fraction of the total particles 
will be condensed to the hubs for 6 < S c I2ll - l23ll . Generally, 
the particle density p = N/L. In the mean field approach, 
the description of for each node is replaced with the mean 
occupation number nrik for the nodes with the same degree k, 
i.e., nik is the average of all the rii at the nodes with degree 
k. Hence m& is not necessarily an integer. In this framework, 
it is shown that the number of mean occupation particle in the 
stationary state is IT23ll 



m k — k/k c ,k < k c \ 
m k = (fc/fc c ) 1/5 ,fc>fc c , 



(1) 



where the crossover degree k c denotes the degree for m% = 1, 
and is given by lT35ll 



I l'V~' 5/,5c fnr 8 <r 8 ■ 

1+1/5) 1 h ' max 1 J u <^ u c , 
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Lp(-7+l+l/<5) 

[j} S [lnk max ]^, for 5 = 5 C 



(2) 
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For 5 > S c , there is no condensation, and the crossover degree 
k c is given by 



Ak: 



7 



:(fccT +2 -^ 7+2 ) + 



Ak-i +1 
7 - 1 - 1/6' 



(3) 



To gain further insight into the relaxation dynamics of the 
ZRP systems B35I1 . we first set up the relationship between 
ZRP and traps model. Treating each node as a trap with tem- 
perature T and depth Ek, we obtain a trapping network where 
both T and E k can be determined by the jumping rate p(m). 
As pointed out in our previous work [ 35 ] , the jumping rate has 
different expressions 



p(m k ) 



m k = k/k c 
m 5 k = k/k c 



for 
for 



(4) 



That is, p(mfc) is proportional to the degree k and has the 
same expression for all degree k. As the jumping rate p(m) 
represents the number of particles jumping out of a given node 
per unit time, the rest time r k is m k /p(mk), which is given by 



1, 



Tie 



for k < k c ; 
,1/5-1 



for k > k c 



(5) 



Obviously, larger p(m) corresponds to higher T and smaller 
Eff, and vice versa. In this way, we compare the rest time 
i~k = rnk/p{mk) in ZRP with the trapping time r k = e l3Ek in 
trap model where (3 = 1/T 13711 . It is easy to obtain: For k < 
k c , /3E k = 0; for k > k c , /3E k = ^-log(-^). Letting T = 
i = 6/(1 — 6), then we have 



E k = 0, for k <k c 
E k = log(k) - log(k c ), 



for k > k c - 



(6) 



We see that 6 — > corresponds to T — > 0, the condensa- 
tion threshold 5 C corresponds to T c = 1/(7 — 2), and 6 — > 1 
corresponds to T — > 00. From Eq. (|6]l we see that larger 
k corresponds to larger E k , i. e., deeper trap. In Ref. H37I1 . 
the average rest time t(S) is exactly the average rest time 
(t) = (fc) _1 Y^k ^P(^) T k, indicating that the framework of 
ZRP is equivalent to that of the trap model. In the thermo- 
dynamic limit of L, N — > 00, the average rest time before a 



ho P isr(*) = £-^ 
equation, we have 



dk. Substituting Eq. @ into this 



TiS) = W>C =m * mdk = W <7) 

which depends on 6 through k c . For 6 < 6 C , substituting 
Eq. d2J into Eq. ^} we obtain 



t(S) 



<fc> Lp( — r +l + l/5)J ^ max ' J°' ^ °c, 

-fcljf [In fcmax]S /or 5 = <5 C . 



(8) 



For <5 C < <5 < 1, the crossover degree fc c can be obtained from 
Eq. ©. Especially, for the case of 6 — 1, we have fc c = (k)/p 
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FIG. 1: (color online), (a) The rest time r k vs degree k for 5 = 
0.0, 0.2, 0.5, 1.0. The solid lines represent the slopes s = 1/5 — 1 = 
1.0 and 4.0, respectively, (b) The average rest time t(S) versus the 
temperature T(8). The "squares" represent the numerical simula- 
tions and the "circles" represent the theoretical results according to 
Eq. Q. The inset shows how t(S) changes with S. 



when m kc = 1 M23ll . Substituting it into Eq. (0 we have 
t = 1, which is consistent with the case of random walk. 

To confirm the theoretical results, we make numerical sim- 
ulations in an uncorrected configuration model (UCM) ll38l . 
We first construct an UCM network with size L = 10 4 and 
degree distribution P(k) ~ Ak~ 3 , where A k 13. Other pa- 
rameters of this network are (k) sw 5, k = 3, and k max = 97, 
where ko and k max denote the minimum and maximum de- 
gree of the network, respectively. We set the particle density 
p = 1 and let the particles hop in the network. In the station- 
ary state, we calculate the mean rest time r k at the nodes with 
the same degree k, and figure out their average r. The results 
shown in Fig.Q]are consistent with the theoretical predictions 
in Eqs. <(5j and (0. As shown in Fig. [T](b), the average rest 
time t(6) decreases significantly with the increase of 5 for 
6 < S c , indicating that the diffusion is remarkably slowed 
down in the condensation phase. 

In the stationary state, a normalized particle distribu- 
tion P eq (k), which is defined as the probability for a parti- 
cle to be in any node with degree k, is given by P eq (k) = 
P(k)m k / p, and thus 



Peg(k) 



Ak-< + 1 ' s 



pk,. 



1/8 



for k < k c ; 
for k > k c 



(9) 



On the other hand, the relaxation dynamics in the conden- 
sation phase is an interesting problem ll22l l23ll . Here we in- 
vestigate how P(k,t w ) converges to P eq (k) in the transient 
state. Until now, there is no exact theory for the relaxation dy- 
namics of the ZRP, and almost all the studies have been inves- 
tigated by Monte Carlo simulations. In the transient period, 
the simulated results stir up one co nject ure that the relaxation 
dynamics has a hierarchical nature ||22| l37ll : At first, subnet- 
work with small degree nodes is stable, and then larger degree 
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regions progressively equilibrate. In fact, the nodes with small 
degrees correspond to shallow minima, which take less time 
to explore, while the nodes with large degrees are deep traps 
which take longer time to equilibrate. At time t w , one can sup- 
pose that the nodes with k < k w are "at equilibrium", while 
the nodes with k ^S> k w are still in the random walk regime. It 
turns out that the particle distribution behaves in each regime 
as 



P(k;t w ) 




k ^ ky'f 

ky k ky 
km k. 



(10) 



Like k c in the stationary state, k„ plays the role of the 
crossover degree scale in a subnetwork with the largest de- 
gree k w and the network size L 1 ~ kf^ , that is, 



k v 



kw S ^ 5c for 5 < 5 C ; 



(In k u 



for 6 = 5 r 



(11) 



Considering that the total time t w is the sum of the trap- 
ping times of the visited nodes, which is dominated by the 
longest one Ty. from Eq. we obtain k v ~ t^"^^^ 1 S ' 
and k w ~ tw for 6 < S c in the looped networks. Fig- 

ure [2] (a) and (b) show indeed that the whole non-equilibrium 
distribution can be cast into the scaling form 



P(k;t w ) = P(k)F (fc/4 M/(1 ~ 5) ) 



(12) 



where F is a scaling function such that F(x) displays the scal- 
ing like Eq. (Q~|i at small x, and F(x) ~ x at large x. In addi- 
tion, as shown in Fig.[2](c) and (d), the whole non-equilibrium 
distribution also displays the scaling 



P(k;t w )=t- s ^G(k/t s w ^) 



(13) 



where G(x) ~ x 1-7 at large x. 

This evolution takes place until the nodes with the largest 
degree k max , equilibrate. For an uncorrected scale-free net- 
work, k max ~ L 1 / 2 so that the equilibration time is 



l.(l-«)/«e 
inax 



L 



1-8 



(14) 



It is consistent with the results in Ref. 11221 12311 . 

The evolution of P(k; t w ) in the condensation phase corre- 
sponds to the aging dynamics of the system, which has a hier- 
archical nature. This dynamics is also investigated through a 
two-time correlation function C(t w + t, t w ) between the states 
of the system at times t w and t w + 1, defined as the the average 
probability that a particle has not changed trap between t w and 
t w + 1 [36]: this amounts to considering that the correlation is 
1 within one trap and between distinct traps. The probability 
that a walker remains in trap i longer than t is simply given by 
exp(— t/ri), so that 



C(tw+t,t w )= / dkP{k-t w )e- t/T ". 



(15) 




FIG. 2: (color online). The scaling form of non-equilibrium distri- 
bution with 5 — 0.0, 0.2 in the transient state. Data in both (a) and 
(b) obey the scaling in Eq. J 121 . those in both (c) and (d) obey the 
scaling in Eq. U3\ . 
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FIG. 3: (color online). C(t w + t,t w ) vs t with S = 0.0 (a) and 
S — 0.2 (b) in the condensation phase. The insets show C(t w + 
t,t w ) vst/t w . 



In Fig. [3] simulations show that the closer ZRP systems ap- 
proach the stationary state, the more slowly particles dif- 
fuse, which seems like the critical slowing down in phase 
transition. As the stationary state is approached, more and 
more particles are trapped in hubs, and it's very difficult for 
these particles to escape from hubs. Thus, the correlation 
length (i. e., the characteristic time for the two-time corre- 
lation to disappear) becomes longer and longer 1I23I . For ZRP 
dynamics in uncorrected scale-free networks, it is obvious 
that the correlation function doesn't obey the simple aging 
C(t w + t,t w ) = g(t/t w ) because of the increasingly parti- 
cle attraction in the transient state. 

Aging properties of the system can also be measured 
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FIG. 4: (color online). The average escape time t esc in both the 
transient state and the stationary state, (a) t eac (t w ) /t%^ c vs t w for 
S = 0.0, 0.2, 0.5, 1.0 in the transient state, (b) t%% c vs 5 where the 
solid line represents the slope s = —logioL — —4.0 from Eq. l !16t . 



through the average escape time t esc (t w ) required by the ran- 
dom walker to escape from the node it occupies at time t w . 
We define t esc = (f) — t w , where t' > t w is the time of 
the first jump performed by the walker after t w , which gives 
tesc{tw) = J dk TkP(k;t w ). For small t w with respect to 
the equilibration time, t esc is growing due to the evolution 
of P(k;t w ). After a long time, mk{t w ) — > in any fi- 
nite system, and then we have t% q c = J/,™"" dk p ^™ kTk , 
Substituting Eq. (Q~|i and (0 into this equation, we can numeri- 
cally calculate t e e q sc for the different 5, t e e q sc 



Jk Q p 

j£ max P J"* 5 — dk. For 6 < S c , the system is in the con- 
densation phase. The nodes with k c < k < k max have 
the capacity to accommodate most particles. Therefore, we 



have^ c 
lows 



P(k) m 2 k - 



-dk, and obtain the scaling as fol- 



t eq ~ k 

b esc ""max 



( 7 -l)(l-5) 



-1-5 



(16) 



In summary, we have set up the relationship between ZRP 
and traps model, i. e., T = 6/(1 — 5), and then investigated the 
relaxation dynamics of ZRP in the framework of traps model. 
The particle interaction, especially the resulted condensation, 
is found to significantly influence the dynamics. In the sta- 
tionary state, a rest time t(5) — pk c /(k) is caused by the 
attraction among particles. In the transient state, a hierarchi- 
cal nature of the aging dynamics for 8 < 5 C is revealed by 
means of the scaling analysis of particle distribution. In ad- 
dition, the equilibration time has the scaling T eq ~ L 1 ^ 5 . 
Moreover, the slow dynamics in both the transient state and 
the stationary state have been compared by two important sta- 
tistical measures, the two-time correlation function and the 
average escape time. Both simulations and scaling analy- 
sis show that the closer ZRP systems approach the stationary 
state, the more slowly particles diffuse. At a long time, the 
average escape time in the condensation phase displays the 
scaling t% q c ~ L}~ s . It is expected that the present work will 
be useful for understanding the slow dynamics of condensa- 
tion in the real world. 
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